arXiv:1501.00106vl [cond-matsoft] 31 Dec 2014 


Effect of Sequence-Dependent Rigidity on Plectoneme Localization in dsDNA 


Shlomi Medalion and Yitzhak Rabin 

Department of Physies and Institute of Nanoteehnology and Advaneed Materials, Bar-Ilan University, Ramat-Gan 52900, Israel 

(Dated: January 5, 2015) 

We use Monte-Carlo simulations to study the effect of variable rigidity on plectoneme formation 
and localization in supercoiled dsDNA. We show that the presence of soft sequences increases the 
number of plectoneme branches and that the edges of the branches tend to be localized at these 
sequences. We propose an experimental approach to test our results in vitro, and discuss the possible 
role played by plectoneme localization in the search process of transcription factors for their targets 
(promoter regions) on the bacterial genome. 


The conformational properties of double-stranded (ds) 
DNA molecules are usually modeled by treating these 
biopolymers as semi-flexible chains with uniform rigidity 
that can be represented by a single persistence length 
(rigidity is proportional to persistence length), Ip 
bOnm in physiological conditions This is justi¬ 

fied when one is interested in large-scale properties of ds¬ 
DNA for which one can replace the sequence-dependent 
distribution of the elastic rigidity by its average over the 
chain. When one is interested in small and intermediate- 
scale phenomena one has to consider the full sequence- 
dependence of the rigidity. While some studies suggested 
that the rigidity of bare dsDNA varies across a limited 
range 40nm < Ip < 7bnm dHHi, experiments on cycliza- 
tion of short dsDNA fragments lOObp) reported much 
higher cyclization ratios than expected (OEj). This led to 
the proposal of DNA kinks - pointlike highly flexible do¬ 
mains © - perhaps due to formation of “DNA bubbles” 
(0®. Such sequence-dependent rigidity is a property 
of bare dsDNA and it has been suggested that the effect 
can be utilized for the design of promoter sequences in 
order to control the DNA binding affinity of transcrip¬ 
tion factors that are sensitive to DNA bendability m- 
It may also arise as the consequence of the binding of 
proteins to specific DNA sequences; indeed, in vivo DNA 
is partially covered by proteins that affect its flexibility 
(|3l [TTHT4]) and/or its local curvature (p!5l [T6|) . For ex¬ 
ample, RecA bacterial proteins polymerize along DNA 
to give an effective persistence length of hundreds of nm 
to the RecA-dsDNA complex ©Ell. Other positively- 
charged proteins (e.g., HMGB) and polyamines, increase 
DNA’s flexibility significantly (pTl ITfi [T8|) . Thus, the 
variability of DNA rigidity may be even higher in vivo 
than that of bare DNA in vitro. 

In this work we study the interplay between local rigid¬ 
ity and plectoneme localization in supercoiled dsDNA. 
When circular DNA is subjected to sufficiently large tor¬ 
sional stress, the minimization of the free energy yields 
strongly-writhed conformations known as plectonemes 
(see e.g. Fig. [^. We show that the number of plec- 
tonemic branches and the locations of the edges (end 
loops of the branches in Fig. are affected by non- 
uniform DNA rigidity. This applies not only to circular 
chains (e.g., bacterial and mitochondrial DNA, and DNA 



FIG. 1. (color online) Representative conformation of a su¬ 
percoiled chain of length of 54005p with a pair of low and high 
rigidity values {lp\l^^) — (55nm, 75n?n). Here Ip^ occupies 
1/3 of the chain and is divided into 4 equally-spaced domains 
of identical size, l\^^ occupies the remaining 2/3 of the chain, 
and the superhelical density is cr = —0.05. The red segments 
represent the soft domains, while the blue ones correspond to 
the stiffer ones. The end loops of the plectonemes tend to be 
localized at the soft domains. 


plasmids), but to topologically-constrained linear dsDNA 
molecules as well, such as eukaryotic chromosomes that 
are attached to the nuclear lamina. 

We performed Metropolis Monte-Garlo (MG) simula¬ 
tions of topology-conserving worm-like rod (WLR) model 
that accounts for bending and twist elasticity of dsDNA. 
A circular dsDNA molecule was modeled as a closed chain 
of N segments of length 306p ^ lOnm each. The twist 
persistence length of the chain was taken to be Itw = 
lAnm (corresponding to Itw = 7.4 in units of segment 
length). The effective diameter of dsDNA was assumed 
to be 5nm, taking into account the screened electrostatic 
repulsion at physiological conditions (Hanoi)- The degree 
of supercoiling of DNA was characterized by the super¬ 
helix density a that is proportional to the amount of tor¬ 
sional stress per unit length of the molecule doami)- Bac¬ 
terial cells have an enzymatic mechanism that keeps the 
superhelix density of the genome at an almost constant 
value, typically in the range —0.03 < a < —0.09 that 
corresponds to somewhat unwound DNA ([CT |2^ and. 
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FIG. 2. (color online) Histogram of the locations of the centers 
of plectoneme edges (end loops), Pedges^ for a circular chain 
with a — —0.06 and 4 soft domains located around the 23rd, 
68th, 113rd and 158th segments. The rigidities correspond 
to the pairs — (dOnm, 55nm) (blue), (55nm, 75nm) 

(green), and (50nm, 200nm) (red). 

unless stated otherwise, we assumed that a = —0.06. 
For more details about the energy form used for accept¬ 
ing/rejecting the Metropolis MC steps see Sec. SI in the 
SI. 

In the simulations we used the pivot (or crankshaft) 
moves described in detail in ref. ©• The angle of rota¬ 
tion around the pivot was tuned in order to achieve the 
desired acceptance rate of 50%. To make sure that the 
simulation was not stuck in a specific plectonemic con¬ 
formation, we used the linking number inversion method 
described in ref. We ran our simulations for 54006p 

long chains with segments of two different rigidities dis¬ 
tributed along the chain such that most of the chain had 
a higher rigidity (/p^^), and a smaller part of the chain 
had a lower rigidity For each chain conformation 

taken into account, we determined the number and the 
locations of plectoneme edges (end loops). The algorithm 
used for this analysis is described in the SI. The chains in 
the simulations were unknotted and we made sure that 
this topology was fixed during the simulation by calcu¬ 
lating the Alexander polynomials and Vassiliev invariants 
after each move, and rejecting topology-changing moves. 

We first considered the effect of non-uniform rigidity 
distribution on plectoneme localization in the case where 
both the stiff and the soft domains are relatively large (as 
for example in the RecA-DNA complexes). The soft do¬ 
main was chosen to occupy 1/3 of the chain and consisted 
of several domains of identical size, equally-spaced along 
the contour. In order to model sequence-dependent rigid¬ 
ity of bare dsDNA we took {lp \ = (40nm, 55nm) (as 

in ref. 0 ) and (55nm, 75nm) (as in ref. ([3|)). For mod¬ 
eling rigidity induced by protein polymerization along 
DNA we used = (50nm, 200nm) (even though 


polymerization of proteins such as RecA yields rigidities 
in the range of 600nm — 800nm (|3l [14)), the effect of 
Ip^^ ^ Ip^ on the localization of plectoneme edges is al¬ 
ready clear for = 200nm). In Fig. [^we plotted a his¬ 
togram of the locations of the centers of plectoneme edges 
for a chain with 4 soft domains. Even for the pair with 
the smallest difference between rigidities (40nm, 55nm), 
there is noticeable localization of the edges to the softer 
segments of the chain. As the difference between rigidi¬ 
ties increases, the localization becomes much more pro¬ 
nounced. For (lp\l^^) = (50nm, 200nm), the edges are 
almost exclusively located in the soft domains. While 
high superhelix density is a necessary condition for the 
formation of plectonemes m, we found that increasing 
a beyond the plectoneme formation threshold does not 
have a noticeable effect on the localization of the edges 
of these branches (see Fig. S5 in the SI). 

For chains of uniform rigidity the number of plectone¬ 
mic branches was shown to increase with a ©• In order 
to understand how the dependence of the number of chain 
branches on a is affected by non-uniform rigidity, we com¬ 
pared chains with 4 soft domains with uniform chains of a 
single persistence length (weighted average of the pair of 
the corresponding persistence lengths). We observed that 
for the (50nm, 200nm) pair, the presence of soft domains 
leads to two-fold enhancement compared to homogeneous 
chains with an average value of Ip = 150nm (the effect 
is much weaker both for chains with lower rigidities and 
those with a smaller difference between higher and lower 
rigidities- see Fig. S3 in the SI). The origin of this effect is 
that a chain with Ip = IbOnm has to pay a large bending 
energy penalty in order to create a tight loop/edge and, 
as a consequence, the number of edges is small. The pres¬ 
ence of soft domains with Ip = bOnm (even when the stiff¬ 
ness of the rest of the chain is increased), allows the edges 
to be localized in these domains, thus decreasing the 
bending energy penalty, and the number of branches in¬ 
creases significantly. The number of plectoneme branches 
is also affected by the number of soft domains. We calcu¬ 
lated the average number of branches as a function of the 
number of soft domains for (40nm, 55nm), (55nm, 75nm) 
and (50nm, 200nm), where the softer part of the chain 
was divided into 2 — 6 domains of equal length, uniformly 
distributed along the chain (Fig. S4 in the SI). As ex¬ 
pected, only a minor effect is observed in the first two 
cases where the difference between rigidities is small but 
in the case of large rigidity contrast, the number of edges 
increases by nearly 50% as the number of soft domains 
increases from 2 to 6. 

Next, we proceeded to characterize the effect of DNA 
kinks (soft local defects) on the localization and the 
number of plectoneme edges. We simulated 54006p-long 
chains with Ip = 50nm, and with 2 — 6 kinks at which 
the persistence length dropped to a value of 4nm (about 
twice the rigidity of ssDNA, ©ESI)- As can be seen in 
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'^kink 

{P'branch) 

pitot) 

kink 

(nn) 

2 

5.99 

0.332 

1.98 

3 

6.28 

0.399 

2.50 

4 

6.41 

0.441 

2.63 

5 

6.62 

0.510 

3.37 

6 

6.76 

0.544 

3.67 


TABLE L The number of branches, (ribranch), the total prob¬ 
ability of an edge to contain a kink, and the total 

number of edges containing a kink, (nn) = (ribranch) • Pklnk 
as a function of the number of kinks in the chain, Ukink, for 
5400bp long chains (cr = —0.06). 


FIG. 3. (color online) Histogram of the locations of plectone- 
mic edges, Pedges, for chains with 4 kinks located at the 23rd, 
68th, 113rd and 158th dimers, with a = —0.04 (dark blue), 
a = —0.05 (light blue), a = —0.06 (light green), a = —0.07 
(orange) and a = —0.08 (red). 


Fig. i for the 4-kink case (and for 2 — 6 kinks in Fig. 
S6 in the SI), plectoneme edges tend to form at the loca¬ 
tions of kinks. Another effect observed in Fig. [^is that 
in the 4-kink case the fraction of edges that contain kinks 
decreases monotonicaly with |cr|, presumably because of 
the number of branches increases monotonically with \a\. 
In the 6-kink case the decrease in the fraction of kink- 
containing edges is observed only for \a\ = 0.07 — 0.08 
(Fig. S9 in the SI). This concurs with the observation 
(Fig. S8 in the SI) that while in the 6-kink case the 
number of edges is only slightly greater than the number 
of kinks in the chain even for the highest value of |cr|, in 
the 4-kinks case the number of edges exceeds the number 
of kinks already for moderate |cr|s. 

In Table [l| and in Figs. S6-S7 in the SI, we show that 
as the number of kinks, rikink^ increases, the number of 
branches, ribranch^ also grows, but this increment is ac¬ 
companied by an even more rapid increase in the frac¬ 
tion of edges that include a kink, , and in the total 

number of edges containing a kink, nn- Even though 
plectonemes have branches even in the absence of kinks, 
when supercoiled DNA contains kinks, the edges of the 
branches are localized preferentially at the kinks; this re¬ 
duces the bending energy penalty for creating edges and 
promotes the formation of new branches. 

As expected from experience with linear worm-like 
chains, the presence of a kink affects the local fluctuations 
of the chain in which it is imbedded. This effect is clearly 
observed in Fig 4(a) [ where we compare the distribution 
of the end-to-end distance d of a short DNA sequence 
(in a long circular chain) of uniform rigidity, with that 
of a DNA sequence of the same length and rigidity but 
containing a kink at its center. Since the length of the 
sequence is only twice the persistence length of dsDNA, 


both distributions are asymmetric, with a peak closer 
to the fully stretched (lOOnm) state. The introduction 
of the kink shifts the maximum of the distribution from 
80nm to 70nm and results in dramatic enhancement of 
the ends contact probability {d < lOnm), an effect that 
has been invoked to explain the anomalously high cycliza- 
tion probability of short dsDNA molecules ([5]-[7j). The 
difference between the distributions becomes much more 
pronounced in the presence of supercoiling. Since the end 
loops of plectonemes tend to form around the kinks and 
the size of the end loop is of order Ip, we expect the ends 
of the lOOnm segment to be preferentially localized in 
the stem of the plectoneme and the most probable end- 
to-end separation will be below lOnm. This expectation 
is fully confirmed in Fig. |4(b)| in which we presented the 
corresponding end-to-end distance distributions of super- 
coiled DNA for a = —0.06 (the histograms for a in the 
range between 0 and —0.07 are shown in Fig. SIO in 
the SI). While the distribution for the no-kink case is 
qualitatively similar to that for cr = 0 (Fig 4(a)| ), the 
peak of the distribution of a kink-containing segment is 
shifted towards contact. In principle, this prediction can 
be tested by FRET experiments by attaching donor and 
acceptor molecules to the ends of a short (of order Ip) 
sequence in a long circular dsDNA molecule, provided 
that one is able to design a sequence that has a high 
propensity for forming a kink (|26H2^ . In the absence 
of supercoiling, the donor and acceptor will be separated 
by tens of nanometers and there will be no FRET signal. 
As torsional stress is introduced into the chain, e.g., by 
raising the concentration of intercalators in the solution 
(|3Ql [3T]) . plectonemes will form with end loops localized 
at the kinks and the distance between the donor and 
acceptor will approach the 6nm limit at which energy 
transfer will take place and a FRET signal will be ob¬ 
served. 


This localization effect may also play a significant role 
in the search process executed by transcription factors 
(TFs) that attach to specific sequences in promoter re¬ 
gions of DNA. Since bacterial DNA is supercoiled, DNA 
segments that reside in the stems of the plectonemes 
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(a) 



(b) 

FIG. 4. (Color online) Histograms of distances of the end 
points of a lOOnm long segment in a 54006p plasmid contain¬ 
ing a kink (light blue) and not containing a kink (orange), for 
(a) torsionally relaxed {a = 0) plasmids and (b) torsionally 
stressed {a = —0.06) plasmids. Representative DNA confor¬ 
mations are shown in the corresponding figures (red), with a 
sequence (light blue) containing a kink (sphere). A sequence 
of the same length but without a kink is shown for comparison 
(orange.) 


are pressed against each other and one expects bind¬ 
ing of proteins to these segments to be suppresed. The 
looped ends of plectonemes are more accessible and, com¬ 
bined with the fact that many bacterial DNA-associated 
proteins have a higher binding affinity to bent DNA 
(di Ea [33]), this suggests that TFs may have a higher 
affinity for end loops of plectonemes. Whether nature 
utilizes such effects to facilitate the search of TFs for 
their DNA binding sites depends on whether promoter 
sites tend to be localized at the edges of plectonemes. 
Indeed, many bacterial promoters contain TATA-box se¬ 
quences that are known for their lower rigidity and higher 
tendency to create DNA bubbles (kinks) ([264[29|) and, ac¬ 


cording to our results, such sequences will tend to nucle¬ 
ate plectonemes and to be positioned at their edges. Fi¬ 
nally, we would like to mention that although we are not 
aware of direct experimental proof of our proposed mech¬ 
anism of plectoneme localization to low-rigidity DNA se¬ 
quences, there has been a recent experimental study of 
the effects of torsional stress on stretched linear dsDNA 
in which repeated hopping of plectonemes between spe¬ 
cific locations along DNA was reported. (|34|) . 
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SUPPORTING INFORMATION 


The Energy Form for the MC Process 


In a discrete worm-like rod (WLR) model that accounts for both the bending diversity and twisting elasticity of 
the DNA, the energy of a chain (in units of ksT) with N segments is given by: 


Ewlr 

ksT 


N 

E 

n=l 






(SI) 


where Ip^n is the dimensionless bending persistence length corresponding to the (n, n + 1) dimer of segments, Itw is the 
dimensionless twist persistence length (both measured in units of the segment’s length As), and kn = 2[1 —cos 
is the dimensionless curvature between the nth and the (n + l)th segments. SuJtw,n = ^tw,n—^tw defines the difference 
between the nth twist angle, ujtw.m and its spontaneous value, . The twist angle is the sum of the first and the 
third rotation angles {uotw.n = 7^ + of fho Euler transformation that rotates the nth segment to the (n + l)th 
segment, while the bending angle, ^n,n+i: is the second Euler rotation angle. 

Topologically, supercoiled dsDNA could be described as a closed chain made of two strands infinitesimally closed 
to each other. This kind of a chain obeys two topological constraints, i.e., the closure of each of its strands. Another 
description for these topological constraints is given in terms of the closure of the center-line of the chain, and the 
closure of the cross-sectional plane. These two constraints are fully accounted by the Fuller-White relation (dHSl); 


Lk = Wr + Tw. 


(S2) 


where the writhe LFr is a measure of the deviation of the center-line from planarity, and the twist Tw is proportional 
to the the sum (in the discrete description) over all the twist angles along the chain. The continuum expressions for 
Tw and Wr are given e.g., in (|3]). In their discrete version they take the forms: 


Tw = -y 

27T ^ 


^tw,n 1 


(S3) 


and: 


Wr = Tyw, 

n,m 


n,m') 


(S4) 
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where Wn,m is the discrete pair-of-segments contribution to the VFr, calculated e.g., in (|4)) (Wr depends only on the 
coordinates of the center-line of the chain). 

For a circular chain, Lk is a topological invariant (an integer) that does not change as long as the topology of the 
chain is maintained. A supercoiled dsDNA is overwound or unwound depending whether its Lk is higher or lower 
than the preferred value of overall twist, = N/27r. When the value of this difference, 

ALk = Lk — differs from zero, the chain becomes torsionally stressed and attempts to minimize its total free 

energy by transforming some of its twist energy into bending energy and, when a critical value of torsional stress is 
reached, plectonemes appear in the chain. 

Since the twist angles of different segments are independent, one may average the twist along the chain by calculating 
{^^tw,n) — 27 tTwIN = 27r{Lk — Wr)/N, where Lk is a pre-determined topological constant, and the Wr is calculated 
using the center-line coordinates of the chain. The resulting twist energy, using ATw = Tw — = ALk — Wr^ 

is then given by: 


Et — 


27t‘^1, 


tw 


N 


ATw^ = 


2nH 


tw 


N 


{ALk - Wrf 


(S5) 


Thus, the twist energy in our simulations is obtained by using only the coordinates of the center-line of the chain for 


calculating Wr, and plugging this value of Wr (and the constant ALk) into Eq. S5 to calculate Et- 


Methods for Identifying Plectoneme Edge Location 

The number and locations of plectoneme edges (end loops) were determined by two methods. In the first analysis 
we used the method of calculation of local writhe, Wr*^^^^^ described in ref Although Wr of a chain is strictly 
defined only for closed chains, one may calculate the Gauss double integral along shorter linear segments of length I 
along the chain to obtain the local writhe associated with the jth monomer, Wr^^^^\j). As shown in Fig. SI, the 
function has maxima at end loops (edges) of the plectoneme, and when the value of such a maximum 

exceeds a threshold, Wr ^threshold) horizontal black dashed line in the figure), this monomer is marked as an edge 
of the plectoneme (a sphere on the chain in the figure). The height of the peaks of Wr^^^^^ depends both on the 
length I of the segments for which we calculate the local writhe and on the superhelix density a, since higher |cr|s 
yield tighter loops resulting in higher local writhe values. Since we want the same threshold value, to 

distinguish between edges to other segments for different as, I must compensate for the a-dependence of the heights 
of the peaks. While in ref the authors flattened the conformations to nearly planar ones before calculating the 
number of edges, we wanted to analyze the three dimensional structures of the plectonemes, and hence, we used a 
different function for /(a) and a different value of In our simulations the choice I = 18/p/(l + 80|a|) 

yielded about the same average heights for the peaks when varying a, and _ q found to be an 

appropriate value to identify the edges from other segments. 

For our second analysis we started from the jth segment, and measured the distance {1) between the (j + l)th 
and (j — l)th segments, the distance (2) between (j + 2)th and (j — 2)th segments, and in general, the distance 
D^{m) between (j + m)th and (j — m)th segments. Plotting D for different js along the chain for a specific realization 
one gets curves such as in Fig. S2. If j is a segment somewhere on the stem of the plectoneme, D will increase almost 
linearly with m as demonstrated, e.g., by the violet curve in Fig. S2, whereas j corresponding to a segment in the 
center of an end loop will behave differently, i.e., it will slightly increase and then will oscillate with a minimal value 
very close to the (excluded) diameter of the chain, d. Running along specific plectonemic conformation (j = 1, 2...), we 
identified the segments that yield local minimum of the sum of distances Ylm=i ^■^ (^) (foi* = 12 and accepted 
only minima under a threshold value of 20 in units of As) as the segments in the center of a loop. An advantage of 
this analysis is that by examining the shape of {m) for the center points of an edge, one may determine the lengths 
of the branches as well, since when the two chain parts depart from each other (for some m), there is an abrupt change 
in the slope of D^{m). However, for our purposes the two methods of analysis agreed in more than 95% of the cases. 
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FIG. SI. (Color online) along the chain, for the specific realization plotted in the hgure. The spheres at the edges 

correspond to the locations of the peaks of 



FIG. S2. (Color online) Typical realization of (m) starting from different js along the chain: center of an edge (orange), 
middle of a branch body (violet), and a point somewhere in between (light blue). 


(5) A. V. Vologodskii, S. D. Levene, K. V. Klenin, M. Frank-Kamenetskii, and N. R. Cozzarelli, Journal of molecular biology 
227, 1224 (1992). 
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FIG. S3, (color online) Average number of branches of the plectonemes as a function of the absolute value of cr, for chains with 
high/low rigidities compared to chains with a single rigidity given by the weighted average of the matching pair of persistence 
lengths: {Ip) = The pair rigidities and the corresponding weighted averges are: {lp \ — (dOnm, 55nm) (orange 

dashed line with circles) and 50n?n (orange solid line with asterisks); (55n?n, 75nm) (blue dashed line with circles) and 68.3nm 
(blue solid line with asterisks); (50nm, 200nm) (violet dashed line with circles) and 150nm (violet solid line with asterisks). 



FIG. S4. (color online) Average number of plectoneme branches as a function of the number of soft domains for chains with 
persistence-length pairs: (lp\lp^^) — (40nm, 55n?n) in orange, (55n?n, 75nm) in blue, and (50n?n, 200nm) in violet. 
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(a) (b) 

FIG. S5. (Color online) Histograms of the locations of the centers of the plectonemic edges (loops), Pedges, for a circular chain 
with 4 soft regions located around the 23rd, 68th, 113rd and 158th segments for: a = —0.04 (dark blue), a = —0.05 (light 
blue), a = —0.06 (light green), a = —0.07 (orange) and a = —0.08 (red), (a) = (40nm, 55nm), (b) = 

(50nm, 200nm). Note that the probability of localization of the edges to the kinks doesn’t change significantly with a. 
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FIG. S6. (Color online) Histograms of the locations of the centers of the plectoneme edges, Pedges, for a circular chain of 
L ~ 54005p ~ 180se^ with point-like kinks with Ip = 4nm where the rest of the chain has Ip = 50nm. (a) two kinks located at 
j = 45 and j = 135 (b) three kinks located at j = 31, j = 91 and j = 151 (c) four kinks located at j = 23, j = 68, j = 113 and 
j = 158 (d) five kinks located at j = 19, j = 55, j = 91, j = 127 and j = 163 (e) six kinks located at j = 16, j = 46, j = 76, 
j — 106, j = 136 and j = 166 (f) the overall probability for an edge to contain a kink as a function of the number of kinks. 
The existence of higher number of kinks effectively “sucks” the probability from the other parts of the chain, and as can be 
seen in (a)-(e) the difference between the probabilities of kink-containing and kink-free edges increases. 



























































Sll 



FIG. S7. (color online) The mean number of plectoneme branches as a function of the number of kinks for 54005p-long chains 
with Ip = 50n?n for the chain and Ip = 4nm for the kinks. 



a 


FIG. S8. (color online) The number of branches of the chain, ribranch^ as a function of |cr| for d-kinks (orange) and 6-kinks (light 
blue). While for 6-kinks the number of branches is comparable with the number of kinks, for 4-kinks the number of branches 
exceeds the number of kinks. 
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FIG. S9. (Golor online) Histograms of the locations of the plectoneme edges, Pedges^ for a = —0.04 (dark blue), a = —0.05 
(light blue), a = —0.06 (light green), a = —0.07 (orange) and a = —0.08 (red), (a) 4 kinks located at the 22nd, 62nd, 112nd 
and 152nd dimers and (b) 6 kinks located at j = 16, j = 46, j = 76, j = 106, j = 136 and j = 166. 
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(a) (b) 




(c) (d) 




(e) (f) 




(g) (h) 

FIG. SIO. (Color online) Histograms of distances of the end points of a lOOnm long segment in a 5400bp plasmid containing a 
kink (light blue) and not containing a kink (orange), for (a) a = 0 (b) cr = —0.01 (c) a = —0.02 (d) a = —0.03 (e) cr = —0.04 
(f) a = —0.05 (g) a = —0.06 (h) a = —0.07. The shift of the maximal probability to short distances, and the significant 
enhancement of probability for the shortest distances are clearly observed for segments containing a kink. 


































































































































































































































































